o 



o 



(N 



in 



BACKSTEPPING DESIGN FOR INCREMENTAL STABILITY 



MAJID ZAMANI AND PAULO TABUADA 



Abstract. Stability is arguably one of the core concepts upon which our understanding of dynamical and 

control systems has been built. The related notion of incremental stability, however, has received much less 

attention until recently, when it was successfully used as a tool for the analysis and design of intrinsic observers, 

^ ' output regulation of nonlinear systems, frequency estimators, synchronization of coupled identical dynamical 

Q ' systems, symbolic models for nonlinear control systems, and bio-molecular systems. However, most of the 

h-^ ^ existing controller design techniques provide controllers enforcing stability rather than incremental stability. 

^H , Hence, there is a growing need to extend existing methods or develop new ones for the purpose of designing 

f**^ ' incrementally stabilizing controllers. In this paper, we develop a backstepping design approach for incremental 

stability. The effectiveness of the proposed method is illustrated by synthesizing a controller rendering a 

synchronous generator incrementally stable. 



1. Introduction 



Stability is a property of dynamical systems comparing trajectories with an equilibrium point or with a partic- 
ular trajectory. Incremental stability is a stronger property comparing arbitrary trajectories with themselves, 
rather than with an equilibrium point or with a particular trajectory. It is well-known that for linear sys- 

^ , tems incremental stability is equivalent to stability. For nonlinear systems, incremental stability is a stronger 

[--». ' property requiring separate concepts and techniques for its study. 



The notion of incremental stability has a long history that can be traced back to the work of Zames in the 60's, 



^+ ■ as described in [Zam96| . In |Zam63| . incremental stability is introduced and studied under the input-output 

setting where control systems are regarded as operators mapping input signals to output signals. Incremental 
f— ^ ■ stability then arises naturally by considering the Lipschitz constant of the operator. A modern treatment 

f^ I of incremental stability, based on Lyapunov methods, appeared only recently in |Ang02| where incremental 

global asymptotic stability ((5-GAS) and incremental input-to-state stability (^-ISS) were defined in a state- 
space setting. The notion of (5-GAS was defined by requiring the Euclidean distance between two arbitrary 
. , - system trajectories to convergqj to zero. Lyapunov characterizations of (5-GAS and S-ISS were also given in 

r> I |Ang02| . There are two other stability properties related to incremental stability that have an equally long, if 

j^ ■ not longer, history. 

The first is the notion of convergent system which, according to |PPvdWN04j . was introduced in the 60 's 
by B. P. Demidovich in |Dem61[ IDem67| . A system is convergent if all the system trajectories converge^ 
to a bounded trajectory. Furthermore, Demidovich also introduced a sufficient condition for a system to 
be convergent, called the Demidovich's condition in |PvdWN05) . Since incremental stability requires every 
trajectory to converge to every other trajectory, an incrementally stable system is also a convergent system 
whenever a bounded trajectory exists. 

The second stability property is contractivity and was introduced in the control community by Lohmiller and 
Slotine in [LS98 although it had been studied before in the mathematical community |Jou05] . Rather than 
comparing trajectories, the notion of contracting system is infinitesimal and requires the decrease of a suitable 
quantity, defined through a Riemannian metric, along trajectories. The definition of contracting system can 
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be seen as a generalization of the Deniidovich's condition. Note that while the Demidovich's condition was 
introduced as a sufficient condition for a system to be convergent, its generalization in |LS98) was directly 
used as the definition of contracting system. 

Although both contractivity as well as convergence are coordinate independent properties, this is no longer 
the case with incremental stability. In this paper, however, we work with a variation of incremental stability 
that is coordinate invariant. This is achieved by no longer insisting on the distance between trajectories being 
measured by the Euclidean metric. 

The number of applications of incremental stability has increased in the past years. Examples include intrinsic 
observer design |AR03) . consensus problems in complex networks |WS05j . output regulation of nonlinear 
systems [PvdWNOS] , design of frequency estimators [SK08b , synchronization of coupled identical dynamical 
systems |RdBS09] , construction of symbolic models for nonlinear control systems [PGT081 IPTOQi IGPT09| , and 
the analysis of bio- molecular systems |RdB09) . Our motivation comes from symbolic control where incremental 
stability was identified as a key property enabling the construction of finite abstractions of nonlinear control 
systems [PGTOSl IPTOQI IGPTOQi IZPJTlOj . Hence, there is a growing need for design methods providing 
controllers enforcing incremental stability since most of the existing design methods guarantee stability rather 
than incremental stability. 

Related work includes the recent design results enforcing the convergent system property through the solution 
of linear matrix inequalities (LMIs) [PvdWNOS, PvdWNOTl IvdWPOS] . In contrast, the results presented in 
this paper do not require the solution of LMIs and existence of controllers is always guaranteed. Backstep- 
ping design methods for incremental global asymptotic stabilitjQ for parametric-strict-fccdbacllj systems were 
proposed before in ^JL02. SK08a_, . and ,SK09,. In this paper, we generalize the results in ^JL02, SKOSaj . and 
^K09\ by: 

1) developing a backstepping design method providing controllers enforcing incremental input-to-state 
stabilitjy and not simply incremental global asymptotic stability; 

2) enlarging the class of control systems from parametric-strict-feedback to strict-feedbacly form. 

The proposed approach was inspired by the original backstepping method described, for example, in [KKK95] . 
Like the original backstepping method, which provides a recursive way of constructing controllers as well as 
Lyapunov functions, the approach proposed in this paper provides a recursive way of constructing controllers 
as well as contraction metrics. Our design approach is illustrated by designing a controller rendering a 
synchronous generator incrementally stable. 

2. Control Systems and Stability Notions 

2.1. Notation. The symbols M, M+ and Wq denote the set of real, positive, and nonnegative real num- 
bers, respectively. The symbols Im, and Om denote the identity and zero matrices in R™x™. Given a 
vector X G K", we denote by Xi the i—th element of x, and by ||a;|| the Euclidean norm of x; we recall 
that ||x|| — \J x\ + a;! + ... + x^. Given a measurable function / : M(j" -^ M", the (essential) supremum of / 
is denoted by ||/||oo; we recall that ||/||oo := (ess)sup{||/(t)||,i > 0}. A continuous function 7 : Rj -^ Rj, is 
said to belong to class K, if it is strictly increasing and 7(0) = 0; 7 is said to belong to class /Coo if 7 £ /C 
and 7(r) ^> 00 as r ^> 00. A continuous function /3 : M(j" x W^ — > K^j" is said to belong to class KL if, for 
each fixed s, the map /3(r, s) belongs to class /Coo with respect to r and, for each fixed nonzero r, the map 
/3(r, s) is decreasing with respect to s and /3(r, s) ^^ as s ^' cxd. If : R" — >■ R" is a smooth function 
with a smooth inverse, called a diffeomorphism, and if X : M" — > R" is a smooth map, we denote by i^^X 
the map defined by ((/)*X)(y) = f^|^ .-w )^ ° <\>^^{v)- Let now G : R" ^ R"^" be a smooth map. The 



"Understood in the sense of Definition l2.2l 
See equation l|3.1|l or [KKK951 for a definition. 
Understood in the sense of Definition [23] 

°See equation I I3.18I I or |KKK95l for a definition. 
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notation (t)*G : M" -^ R"^" denotes the smooth map {(f>*G){x) = {^VG{(f>{x)){^). A Riemannian metric 
G : K" -^ R"^" is a smooth map on R" such that, for any x G W\ G{x) is a symmetric positive definite 
matrix [Lee03j . For any x G M" and smooth functions /, J : M" -^ R", one can define the scalar function 
(7, J)g as l'^{x)G{x)J{x). We will stiU use the notation (/, J)g to denote I^GJ even if G does not repre- 
sent any Riemannian metric. A function d : R" x R" — > R^j" is a metric on R" if for any x,y,z G R", the 
following three conditions are satisfied: i) d{x,y) = if and only ii x — y; ii) d{x,y) = d{y,x); and iii) 
d{x, z) < d{x, y) + d{y, z). We use the pair (R", d) to denote a metric space R" equipped with the metric d. 
We use the notation dg to denote the Riemannian distance function provided by the Riemannian metric G 
|Lee03) . 

2.2. Control Systems. The class of control systems that we consider in this paper is formalized in the 
following definition. 

Definition 2.1. A control system is a quadruple: 

I] = (R",U,W,/), 
where: 

• R" is the state space; 

• U C R™ is the input set which is compact and convex; 

• U is the set of all measurable functions of time from intervals of the form ]a, 6[C R fo U with a < 
and b > 0; 

• / : R" X U — !> R" is a continuous map satisfying the following Lipschitz assumption: for every compact 
set Q C R", there exists a constant Z G R^ such that \\f{x, u) — /(y, u)|| < Z\\x — y\\ for all x,y ^ Q 
and all u ^ U. 

A curve ^ :]a, b\-^ R" is said to be a trajectory of E if there exists v €U satisfying: 

(2.1) m^f{at),v{t)), 

for almost all t G ]a, 6[. We also write Cxv{t) to denote the point reached at time t under the input v from 
initial condition x = ^x-u(O); this point is uniquely determined, since the assumptions on / ensure existence 
and uniqueness of trajectories |Son98j . A control system S is said to be forward complete if every trajectory 
is defined on an interval of the form ]a, oo[. Sufficient and necessary conditions for a system to be forward 
complete can be found in |AS99) . A control system E is said to be smooth if / is an infinitely differentiable 
function of its arguments. 

2.3. Stability notions. We start by introducing the following definitions which were inspired by the notions 
of incremental global asymptotic stability ((5-GAS) and incremental input-to-state stability ((5-ISS) presented 
in |Ang02| . 

Definition 2.2. A control system E is incrementally globally asymptotically stable (S^-GAS) if it is forward 
complete and there exist a metric d and a ICC function f3 such that for any t G R^ , any x,x' G R" and any 
V E U the following condition is satisfied: 

(2.2) d{i,^{t),^,,^{t)) < p {d{x,x') ,t) . 

While 5-GAS, as defined in |Ang02| , requires the metric d to be the Euclidean metric. Definition 12.21 onlv 
requires the existence of a metric; hence, the existential quantifier in the acronym (5g-GAS. We note that while 
(5-GAS is not invariant under changes of coordinates, (Jg-GAS is. If (p : R" -^ R" is a bijective change of 
coordinates, inequality (|2.2|) transforms under to: 

d' (0 o Cx^W, </> ° ^.'v{t)) < /3 (d' (0(x), </)(x')) , t) , 

where d'{y,y') = d{<j)~^{y),<j)~^{y')). We shall return to the comparison between Jg-GAS and 5-GAS in 
Subsection l2.4l Nevertheless, when the origin is an equilibrium point for E and the map ip : R" -^ Rj, defined 
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by i^{x) — d(a;, 0), is continuous and radially unboundeqj, both (J^-GAS and (5-GAS imply global asymptotic 
stability. 

Definition 2.3. A control system S is incrementally input-to- state stable (8^-ISS) if it is forward complete 
and there exist a metric d, a ICC function (3, and a ICoo function 7 such that for any t € M.q , any x, x' £ M", 
and any v, v' € U the following condition is satisfied: 

(2.3) d (6. W, 6'.' (t)) < /? (d {x, x') , i) + 7 (lit; - v'WJ . 

By observing p.2p and (|2.3|) . it is readily seen that ^a-ISS implies ^a-GAS while the converse is not true in 
general. Moreover, whenever the metric d is the Euclidean metric, i5a-ISS becomes (5-ISS as defined in |Ang02| . 
We note that while 5-lSS is not invariant under changes of coordinates, (5g-ISS is. Once again, although ^g-ISS 
is not equivalent to 5-ISS, both notions imply input-to-state stability whenever the origin is an equilibrium 
point for S and the map ?A : R" — > Mq , defined by '4!{x) — d(a;, 0), is continuous and radially unbounded. 

2.4. Descriptions of incremental stability. One of the methods for checking incremental stability proper- 
ties consists in using Lyapunov functions. The Lyapunov characterizations of 5-GAS and (5-ISS were developed 
in |Ang02| . In this paper we follow an alternative approach based on contraction metrics describing (5g-GAS 
and Jg-ISS properties, when d is the Riemannian distance function. The notion of contraction metric was 
popularized in control theory by the work of Lohmiller and Slotine JLS98J that relies on variational systems. 

The variational system associated with a smooth control system S = (M", U,W, /), when we only have state 
variations, is given by the differential equation: 



(-) >^S 



5L 



for any v £ U and where SS, is the variatiorO of a trajectory ^ of S. More details about the variational 
system can be found in |CS87| . Similarly, the variational system associated with a smooth control system 
S = (M", U, W, /), when we have state and input variations, is given by the differential equation: 

(2.5) ^M)='/ ^e + ^ Sv, 

at ox x=5 ou ^=e 

where S^ and 6v are variations of a state trajectory f and an input curve u of E, respectively. 

The following definition is adapted from |LS98I : 

Definition 2.4. Let S = (R", U,U, f) be a smooth control system equipped with a Riemannian metric G. The 
metric G is said to be a contraction metric, with respect to states, for system S if there exists some A G M"*" 
such that: 

(2.6) {X,X)f<-\{X,X)g 

for F{x,u) = (^\ G{x) + G{x)^ + ^f{x,u), anyueU,Xe M" and x G M", or equivalently: 

(2.7) ^^((^) G{x) +G{x)^ + ^fix,u)\ X < -XX^G{x)X, 

where the constant A is called contraction rate. 

When the metric G is constant, the condition (|2.6I) or (|2.7p is known as the Demidovich's condition [PvdWNOS] . 
It is shown in |PvdWN05) that such condition implies incremental stability and the convergent system property. 



Under the stated assumptions it can be shown that a(||a;||) < ^(x) < a(||a;||) for /Coo functions a and a. 
'The variation S^ can be formally defined by considering a family of trajectories S^xv (t, e.) parametrized by e G R. The variation 



of the state is then 5^ = g " . 
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Note that the inequality (|2.6p or (|2.7p imphes: 

(2.8) j^{S^,SOg<-X{S^,SOg, 

when we only have state variations and S^ is the variation of a state trajectory ^ of S. 
The following theorem shows that the inequality (|2.6p implies (5g-GAS. 



Theorem 2.5. Let E = (R", [),U,f) be a smooth control system equipped with a Riemannian m,etric G. If G 
is a contraction metric, with respect to states, for S and if (M",dG) is a complete metric spaced, then S is 
dj-GAS. 

Different variations of this result appeared in |LS98j and JAR03) ; see |AR03| for a concise proof and |ZPJT10] 
for a more detailed treatment including the completeness assumption. It is also shown in }ZPJT10] that the 
additional assumption ui{X, X)j^ < {X, X)a < ^{X, X)i^ for cj, oJ G M^ leads to the stronger conclusion that 
E is in fact S-GAS. 

If, in addition to state variations, we also allow for input variations we obtain the notion of contraction metric 
with respect to states and inputs. 

Definition 2.6. Let E = (R", \J,U,f) be a smooth control system on R" equipped with a Riemannian metric 
G. The metric G is said to be a contraction metric, with respect to states and inputs, for system E if there 
exists some A € R^ and a e R^j" such that: 



(2.9) {X,X)f + 2(^Y,X) <-\{X,X)G + a{X,X)l{Y,Y)l 



\d 



G 



for F{x,u) = ffl^ G{x) + G{x)^ + ^f{x,u), any X e R", a; G R", u £ U, and Y G M™, or equivalently: 
X^ U^YG{x)+G{x)f^ + ^f{x,u)) X + 2Y^ {I{:Yg{x)X 



(2.10) < -\XTG{x)X + a{X^Gix)X)HYTY)^ , 
where the constant A is called contraction rate. 

Note that the inequality ([2^ or (|2.10p implies: 

(2.11) 1^'^^''^^^« - -MS^,S0G + a{S^,S0i{Sv,6v)]^, 

when we have state and input variations and where J^ and Sv are variations of a state trajectory ^ and an 
input curve u of E. 

The following theorem shows that the inequality (|2.9p implies (5g-ISS. 

Theorem 2.7 ([ZPJTTO]). Let E = (R", U, W, /) be a smooth control system on R" equipped with a Riemann- 
ian metric G. If the metric G is a contraction metric, with respect to states and inputs, for system E and 
(R",d(3) is a complete metric space, then E is S^-ISS. 

Similarly to contraction metrics with respect to states, it shown in |ZPJT10] that E is S-ISS when the inequal- 
ities a;(X, X)/„ < {X,X)g < uj(X,X)i^ are satisfied for w, ZJ G R+. 

In the next section, we propose a backstepping design procedure to render control systems incrementally 
stable. 



A metric space (M", d) is said to be complete if every Cauchy sequence of points in R" has a limit that is also in K". 
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3. Backstepping Design Procedure 

The method described here was inspired by the original backstepping described, for example, in IKKK95J . 
Consider the class of control systems E = (M" , U , W, /) with / of the parametric-strict-feedback form KKK95 : 

fi{x,u) = hi{xi)+biX2, 

f2{x,u) = h2{xi,X2) +b2X3, 

(3.1) : 

fn-l{x,u) ^ hn^i{xi,--- ,Xn-l) + bn-lXn, 

fnix,u) = h„{x)+g{x)u, 

where x E K" is the state and u G U C M is the control input. The functions /i^ : R* — > R, for i = 1, . . . , n, 
and g : R" — ?► R are smooth, g{x) ^ over the domain of interest, and hi E R, ior i = 1, . . . ,n, are nonzero 
constants. 

We can now state one of the main results, describing a backstepping controller for control system (13. ip . 

Theorem 3.1. For any control system E — (R", (J,U,f) with f of the form iS. 1]) and for any A £ R+, the 
state feedback control law: 



(3.2) 
where 



k{x,u) 



9{x) 



kn{x) - hn{x) 



1 



9{x) 



ki{x,u) = -bi-i {xi-i - (l)i^2{x)) - -z{xi- (j)i^i{x)) -\ -^f{x,k{x,u)), for I = !,■■■ ,n, 



Mx) 



, for 1^1,- ■■ ,n-l, 



— ki{x) - hi{x) 

(t>-i{x) = 0o(a;) = Vx e M", 6o = 0, and xq = 0, 
renders the control system E 5^-GAS. 



Proof. Consider the following system: 

where 77; = [^i,--- ,^;_i]^,B; = [0, • • • ,0,&i_i]^ £ R'"\ z; = [yf xi]'^ e M' is the state of Ej, y; = [a;i,--- ,a;;_i]^, 
and Fiijji) = [/i(x, u), • • • , fi-2{x, u), hi^i{xi, • • • , xi-i)]'^ . By using induction on I, we show that the metric 
G/, defined by: 



(3.4) 



Giiyi) = 






is a contraction metric, with respect to states, for the system p.3p with contraction rate A. For / = 1, it can 
be easily checked that Gi(yi) = 1 is a contraction metric, with respect to states, with the contraction rate A 
for the scalar system: 

Ei:a = fci(a) = -^a- 

Assume that the metric Gk~i is a contraction metric, with respect to states, for the system Efc_i, for some 
2 < k < n, and with contraction rate A. This implies: 

T 

V ovk J ayu 



y- x\ ( r^^(^k±M^V c,_,(.,_,) + G,_,(.,„,)^^(^:k±M^ -f ^^ (F. + B,,,.,) 



(3.5) 



dVk 
<-\\y'^ x\Gk-i{yk-i) 



Y 
X 



Y 
X 
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dzk 



(3.9) [Y^ X] { I "^^'-^""^t^"'""' ) G.(y.) + G.(y.) "'"'"+^%1'"°'"°']" + G.iy,, 



Y 
X 



Y' X] 



(Fk+BkXk) 



9£fc 
dVk 



Syi"^ + 



d<t>k 



( 8tt>k-i\ 84>k-i 

SVk 



8<l>k 



dVk Svk 2 dy 

T 



Bk 

- ~ BiGk-iivk-i) ~2+ a«7^ -^fc 



9»fc 



8!/t 



1 

9£fe 
dVk 



+ 



G.-i(y.-i)+(^) 



S<*fc- 



(Ffc + BfcXfc; 



T S^<iifc_i S0fc^i SFfc , A 9<t>k-l 



dy. 



"2 H 



Syi- Safe 2 Si/i. 



SSfc 



BlGk-i(yk-i) 



9yk 



Bk 

'2 + Syfc 



ififc 



1 

+ 



""^ [Fk + B.xO + %^(F. + B.x.)%l + i'-^YiFk + BkXkf'-^ -^^(F. + BkXk) 



9yk 



' 3(F^+Bk^k~l) 
dyk 



-{Fk + BkXk 



9yk \ 8yk 

T8^<t>k~ 



'9w7 

Y ' 
X 



1 Gk-i{yk-i) + Gk~i{yk-i) 



a(Fk + Bk4>k-l) , 8Gk-i 



8vk 



+ 



8yk 



Fk]-\ 



,-i\'^ d4>k-i ^/84'k-i^'^ 



8vk 



Y 
X 



8vk 

< -X [F^ X] Gkivk) 



8yk \ dyk 

-X 



Y 

X 



for any F £ M , and X e M. Since the metric Gk-i is only a function of yk-i = [xi,- ■ ■ , 2:^-2] , and the 
vector Bk has zero entries except for the last entry, it can be easily shown that g''~^ Bk — Ok-i, and the 
inequality p.Sp reduces to: 

'r d{Fk + Bk(l^k-i '^'^ 



dyk 



(3.6) 
Now, we show that: 

(3.7) 



r ( \^r I. , diFk + Bk<t>k-i) , OGk-i ^ 

Gk-i{yk-i) + Gk-i(yk-i) 7{ 1 — 7^ -Tfc 

oyk oyk 

< -A [r^ x] Gk^yk-i 



Y 
X 

Y ■ 

X 



Gkivk) 



G 



k-i{yk-i) + [-^^ j -^^ -[^^) 

dct>k-i 1 



dVk 



is a contraction metric, with respect to states, for the system S^. For any nonzero vector [y-^ X\ G R'^, we 
have: 



[y^ X] Gk{vk) 



Y 
X 



[Y-X] 



G._i(,._,) + (^)' 



dVk 



\ dyk J 



dyk 



Y 
X 



(3.8) 



Y^Gk-i{yk-i)Y + 



d<j)k- 
dyk 



-Y -X 



If y e M'^^^ is the zero vector, X must be nonzero implying that the equation p.Sp is equal to X'^ which is 
positive. On the other hand, if y S R*-'"^ is nonzero, Y'^Gk-i{yk-i)y is a positive scalar because Gk-i is a 
Riemannian metric. Hence, Gk is positive definite. Using the inequality p.6p . the long algebraic manipulations 
in p.9p show that Gk satisfies (|2.7p with the contraction rate A. Hence, the metric Gk is a contraction metric, 
with respect to states, for the system E^. Therefore, for any / < n, the metric Gi is a contraction metric, with 
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(3.10) 



respect to states, for the system (13.31) and with the contractfon rate A. 

The proposed control law (|3.2p . transforms a control system of the form (|3.1I) into: 

^n = kn(Tjn, £,n) + V. 

It can be easily checked that v does not appear in the variation of I]„ when we only have state variations. 
Since the metric Gn is not a function of the n-th state, its derivative with respect to time does not include v. 
Hence, we can apply the induction results to I]„ to conclude that the metric G„ is a contraction metric, with 
respect to states, for I]„ and with the contraction rate A. Moreover, it can be readily seen that Gn = 4'* In, 
where 

Xi 

X2 -(t)i{x) 
X3 -02 (a;) 



(3.11) 



ip{x) = 



Xn - '/'n-l(a;) 

Note that d/^ is just the Euclidean metric and we know that (R", d/^) is a complete metric space. Moreover, 
since ip • J^" ~^ H^" is an isometrjij, (R", d^^) is also a complete metric space [Lee03' . By using Theorem l2.5l 
we conclude that a control system of the form p.ip . equipped with the state feedback control law p.2p . is 
Jg-GAS. The ^g-GAS condition ((2?2|) . as shown in }AR03| . is given by: 



(3.12) 



{(xv{t),^^'v{t)) <e~^'dG„ix,x'). 



a 



Remark 3.2. The contraction metric G„ (?/„), with respect to states, for the control system iS.l]) . equipped 
with the state feedback control law i3.2\) . is given by: 
(3.13) 







- 


[11 1 f S'l>i\ S<t>i 
' ' '^ \9y2J ay2 

a<l>^ 
L ai/2 


(a<t, 

\dy 

1 

302 

8y3 













d4>2 \ 302 
^ 9y3 J 9y3 



9<Pn-l 



9<t>2 
, SH3 , 



' d4>„- 



a<t>„ 



' 9rt>„- 



dy-n 



6yn 



dyn 



Sy„ 



where yi — [xi,- ■ ■ , xi-i]"^ , for I — 2, ■ ■ ■ ,n, and the contraction rate is A. 
Remark 3.3. It can be checked that the function 

^(^) = 2 ^^^'+^ -Mx)f , 
1=0 

is a Lyapunov function |Kha96j for the control system i3.1\) . equipped with the state feedback control law iS.^) 
when u = 0. Moreover, the Hessian of V{x) is equal to the contraction metric Gn, with respect to states, 
defined in i3.13\) . 



In the next theorem, we show that control law 

Theorem 3.4. For any control system T, — (l 
state feedback control law: 



(3.14) 



k{x,u) 



9{x) 



also enforces (5g-ISS. 
\ U, W, /) with f of the form fiS. 1]) and for any A G R+, the 



kn{x) - hn{x) 



9{x) 



Suppose M and A/ are Riemannian manifolds with Riemannian metrics G and G, respectively. A smooth map »/> : M — > M 
is called an isometry if it is a diffeomorphism satisfying G = tp'O. 
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where 



ki{x,u) 



-bi^i {xi^i 



(f>i-2{x)) - - {xi -0/_i(x)) + g^ f{x,k{x,u)), for I 



, for 1^1,- ■■ ,n- 1, 



,n, 



— ki{x) - hi{x) 

(t)-i{x) = 0o(x) = Vx e M", bo = 0, and Xq = 0, 
renders the control system E S^-ISS with respect to the input v. 



Proof. Consider the following system: 

(3.15) S, : f,i^ Fi{rii) + Bi(j>i-i{r]i), for / = 2, • 



,n, 

iT 



where rji = [Ci,--- ,6-i] : S; = [0, • • • ,0,6i_i]^ e R' \ yf = [a^i,--- ,xi^iY' is the state of E;, and Fi{r]i) 
[/iCCj"*^): ■ ■ ■ 1 /i-2(Ci ^)5 ^i-i(Ci: ■ ■ ■ T^i-i)]'^- As proved in Theorem lOl the metric Gi in p. 41) is a contraction 
metric, with respect to states, for the system E; and with the contraction rate A. The proposed control law 
p.l4p . transforms a control system of the form (|3.ip into: 



(3.16) 

Now, we show that: 

(3.17) 






GniVn) 



G. 



n-l{yn-l)+ [^^) 



dyn 



\ dy„ ) 



dVn 



is a contraction metric, with respect to states and inputs, for the control system (j3.16p . For n = 1, it can be 
easily checked that Gi{yi) = 1 is a contraction metric, with respect to states and inputs, with the contraction 
rate A, satisfying (|2.10p with a = 2 for the scalar control system: 



s:a = fci(a) 



^ = -2^1 



As proved in Theorem 13. 1[ GniVn) is positive definite. Using the inequality (|3.6p for k = n, long algebraic 
manipulations similar to those in p.9p show that G„ satisfies (|2.10p with the contraction rate A and a = 2. 
Hence, the metric Gn is a contraction metric, with respect to states and inputs, for the control system p.l6p . 
As explained in the proof of Theorem 13.11 we know that (]R",dG„) is a complete metric space. By using 
Theorem 12.71 we conclude that a control system of the form p.ip . equipped with the state feedback control 
law p.l4p . is J^-ISS with respect to v. The Jg-ISS condition (12.31) . as shown in [ZPJTIO^ . is given by: 

dGA^xv{t),^x'v'{t)) < e-**dG„(x,x') + |(l-e-**) |10-9'||oo 



< e 2*dG„(a;,x') + yI|w - w 



Remark 3.5. The contraction metric, with respect to states and inputs, for a control system of the form 
113.1]) . equipped with the state feedback control law ^3.14^ , is given by 113. 13\) . 

Remark 3.6. It can be shown that the function 



-. n— 1 



1=0 



is an Input-to- State Stability Lyapunov function |Kha96| with respect to v for a control system of the form 
113. 1\) . equipped with the state feedback control law |j'. J^p. 
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Now, we extend the results in Theorems 13.11 and 13.41 to the class of control systems E — (M",R, W, /) with / 
of the strict- feedback form |KKK95| : 

fi{x,u) = hi{xi) + gi{xi)x2, 

f2{x,u) = h2ixi,X2) + g2{xi,X2)x3, 

(3.18) : 

fn-l{x,u) = hn~i{xi,--- ,Xn-l) + gn-lixi,- ■■ ,Xn-l)Xn, 
fn{x,u) = hn{x) + gn{x)u, 

where x G M" is the state and m S U C M is the control input. The functions /i^ : R' — ^ K, and gi : M* — !■ R, 
for i = 1, . . . ,n, are smooth, and gi{xi, • ■ • , Xi) ^ over the domain of interest. 

In order to extend Theorems 13. 1[ and l3.4l to control systems of the form (|3.18p . we need the following technical 
lemmas. 

Lemma 3.7. Let E = (M", U,iY, /) he a control system and let cj) : R" — ?> R" be a smooth map with a smooth 
inverse. If the metric G is a contraction metric, with respect to states, for T,' = (R", U,W, (/>,/) and with 
contraction rate X £ R+, then the metric 4>*G is a contraction metric, with respect to states, for the system E 
and with the contraction rate A. 

Proof. Since G is a contraction metric, with respect to states, for the system E' and with the contraction rate 
A, using the inequality (|2.8p . we have: 



(3.19) j^{S7^,Sv)G<-HSr,,Sv)G, 

where Srj is variation of the state trajectory of E'. Since G is a metric and Q{x) — -if-{x) is an invertible 
niatri?l^°l. it is readily seen that {(j)*G){x) is a positive definite matrix. We now show that the metric 0*G is a 
contraction metric, with respect to states, for the system E. For the coordinate transformation ry = 0(C): ^^ 
have: 

(3.20) Sij = e(0<5e 

By taking the derivative of (|3.20p with respect to time, we obtain: 

(3.21) ^^s^ = e{OS( + e{Oj^s^. 

Using (13311), (13201), and (P?^ . we obtain: 

— {5ri,5ri)G = (jt^V GSri + 6ifG—5ri + 5ri^G5ri 

Q5S, + Q^si] GQ5S, + {Q5S,f G (q5S, + Q^5d\ + SjfGSr, 
dt J \ dt J 

|<5e) rCS^ + 5e<t>*Gj^5^ + <5e| (r G) 5^ = j^{5^, SO^^-g 

< -X{6t], Sr])G = -XSri'^GST] = -X (QSO'^ GQS^ = -Xd£,'^(l)*GS£, 

= -X{StSO^'G- 

Hence, the metric (J3*G is a, contraction metric, with respect to states, for the system E and with the contraction 
rate A. D 

In the next lemma, we extend the results of Lemma 13.71 to contraction with respect to states and inputs. 



'^For any smooth map (f> : R" — >■ R" with a smooth inverse, it is easy to show that -jt^(x) is an invertible matrix for any 
X GR". 
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Lemma 3.8. Let S = (M", U,W, /) he a control system and let (j) ■ K" -^ K" be a smooth map with a smooth 
inverse. If the metric G is a contraction metric, with respect to states and inputs, satisfying h2.10fj with 
contraction rate A € R+, and a £ Mj for E' ~ (M", U,Z//, </)*/), then the metric (f>*G is a contraction metric, 
with respect to states and inputs, satisfying i2.10\) with the contraction rate X, and the nonnegative constant 
a for the system T,. 

Proof. Since G is a contraction metric, with respect to states and inputs, for the control system E', satisfying 
(|2.10p with the contraction rate A, and a E R^ , using the inequahty (12.111) . we have: 



(3.22) di^^''^^'^^^ - -H6v,Sv)g + a{6v,Sv)l{dv,Sv)l 

Using ((3?20| . ([33T|) . ((3?22]) . and the resuhs of Lemma [3Jl we obtain: 



dt 



{Sri,6r])c 



dt 



{s^,so^'G 



< 



-X{5f],Sr])G + a{Sr],6ri)^{Sv,Sv)f 
Hence, the metric (f>*G is a contraction metric, with respect to states and inputs, satisfying (j2.10p with the 



contraction rate A, and the nonnegative constant a, for S. 



D 



We can now state the main resuh for a control system S = (R", (i,U, f) with / of the form (I3.18p . 
Theorem 3.9. Let E — (K",U,Z//, /) be a control system where f is of the form US. 18]) . The state feedback 



control law u = k{(j){x),u), where k was defined in VS.'2\) and 
inverse) defined by: 



l" is the smooth map (with smooth 



(3.23) (/>(x) 



renders control system E 5^-GAS. 



Xi 

gi{xi)x2 
gi{xi)g2{xi,X2)xz 

YY^Zi 9ii.xi, ■ ■ ■ ,Xi)x 



Proof. The coordinate transformation rj — 0(^) transforms the control system E = (M", U,i^, /) with / of the 
form p.lSp to the control system E' = (R", U, W, /'), where /' — (t)*f . It can be easily checked that /' has the 
following form: 



(3.24) 



f2iy,u) 

fn~i{y^u) 
fniy^u) 



h'liyi) + y2, 
h'2{yi,y2) + 2/3, 

K-i(2/i,--- ,yn-i) 
h'niy) + g' {yh, 



where /i^ : M* — >■ M, for z = 1, ■ • ■ ,n, are smooth functions, g' = YViZi 9i' ^^^ y ^ ^" ^^ ^^^ state of E'. As 
proved in Theorem 13.11 the state feedback control law fc, defined in (|3.2p . makes the metric Gn, defined in 
p.l3p . a contraction metric, with respect to states, for the control system E' and with the contraction rate A. 
As proved in Lemma l3.71 the metric (jfGn is a contraction metric, with respect to states, with the contraction 
rate A, for the control system E, equipped with the state feedback control law k{(f>{x),u). Since (R",dG„) is 
a complete metric space and cj) is an isometry, (M", d^.g^) is also a complete metric space |Lee03| . Therefore, 
the state feedback control law k{(t>{x),u) makes the control system E (5g-GAS. D 
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The ^3-ISS version of Theorem 13.91 is given by the following result. 



Theorem 3.10. Let S — (R", (J,U,f) be a control system where f is of the form f3.18\} . The state feedback 
control law u — k{(j){x),u), where k and (j) were defined in {S.l^-^ and i3.23\) . respectively, renders control 
system S 6j-ISS with respect to the input v. 



Proof. By following the same steps as in the proof of Theorem l3.91 and using Lemma 13.81 we obtain that the 
state feedback control law u — k{(j){x), u) makes the metric (j)*Gn a contraction metric, with respect to states 
and inputs, for S and with the contraction rate A. Hence, the control system S, equipped with the state 
feedback control law k[(j){x), u), is ^g-ISS with respect to the input v. D 



Remark 3.11. Although we only discussed single input control systems, extensions to multi input control 
systems are straightforward using the techniques in JKKK95] . 



4. Backstepping Controller Design for a synchronous generator 



We illustrate the results in this paper on a synchronous generator [RGHSOf] connected through a transmission 
line to an infinite bus. The control system E — {S? , \J,U, f) with / of the form: 



(4.1) 



/i(x,m) = X2, 

f2{x,u) = -Ex2 + FPmQ + VsGeqo sin{So + xi) + VsG sm{do + xi)x3, 
fsix, u) = -1x3 + JVs sin((5o + a;i)a;2 - leqo + IKcU, 



models a synchronous generator connected to an infinite bus. In the aforementioned model, xi is the deviation 
of the power angle, X2 is the relative speed of the rotor of the generator, 3:3 is the deviation of the quadrature 
axis voltage of the generator. So is the operating point of the power angle, Pmo is the operating point of the 
mechanical input power, Cgo is the operating point of the quadrature axis voltage of the generator, Vs is the 
infinite bus voltage, Kc is the gain of the excitation amplifier, and u is the input of the silicon-controlled 
rectifier amplifier of the generator. Other parameters in (14.11) are given by: E 



D 
2H' 



j_ 

T' 



F 



2ff' 



G 



"0 1 

2H X„. 



and J 



x^-x. 



X' 



-, where D is the per-unit damping constant, H is the inertia constant, wq is 



the synchronous generator speed, T' — jt^T^q, Xqs — Xt + hX 



^q^ ^'ds 



Xj 



\Xl + X' T' is the 



direct axis transient short-circuit time constant, Xt is the reactance of the transformer, Xq is the quadrature 
axis reactance, X'^ is the direct axis transient reactance and X^ is the reactance of the transmission line. We 
assume that sin(5o + xi) is nonzero over the domain of the interest. 

The control system (|4.ip is of the form (|3.18p . The coordinate transformation (I3.23p . given by: 



(4.2) 



m 



m 



6 

14Gsin((5o+Ci)e3 



transforms the control system S = {S? , \J,U,f^ to the control system S' = 
the form: 



— fia3 



\J,U,f) with /' = <^,/ of 



f[{y,u) = h[{yi) + y2 ^ y2, 
(4.3) f2{y,u) = h'2{yi,y2)+y3 = -Ey2 + FP^o + VsGeqosm{So+yi)+y3, 

fsiy^u) = h'3{y)+g'{y)u = -IVsGeqosm{So + yi) + JV,^Gsin^{So + yi)y2 
cot(5o + yi)y2y3 + IKcVsGsm{5o + yi)u. 



Iy3 
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By using the results in Theorem 13.41 for a control system of the form (|4.3p and for A = 2, we have: 

hiv) = {-5 + 3E-E^)Tj2-3m + {E-3)Tj3 + {E-3)FP,nn 
+ {E ~3) VsGeqQ sin(Jo + Vi) - VsGeqQ cos(5o + Vi)^- 
Therefore, the state feedback control law: 

1 ., 1 ^ 
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(4.4) k{r],v) 



.'(,)"=<"'-'■'<"»'.'(,)■ 



1 



IKrVsCsmiSo + Tji) 



-5 + 3E- E^) 772 - 3771 + (S - 3 + /) 773 + (^ - 3) FPmo 



(S - 3 + /) VsGeqQ sin((5o + m) - VsGcgo cos((5o + 771)772 - JV^^Gsin^iSo + 771)7^2 
cot((5o + 771)7^2773 



IK,VsGsm{5o + m)' 



makes the control system S' (5g-ISS with respect to the input v. The corresponding contraction metric, with 
respect to states and inputs, for the control system (|4.3p is given by: 



Giy) 



14_ ( ^^ 
'-^ {.dyi 



901 
8yi 



a<t>i 
9yi 



\dyi , 



902 
dz2 



902 
9Z2 



2 + {2 + VsGeqQCos{6o+yi)f 

-2E + 5-iE- 2) VsGeqo cos(<5o + yi) 

2 + VsGcgo cos(5o + J/i ) 

iT 



-2E + 5 - {E - 2) VsGeqo cos(<5o + yi) 2 + VsGe^o cos(<5o + yi) 
(E- 2)^ + 1 2-E 

2-E 1 



where zj — [7/1 7/2]"^- By using Theorem 13.101 the state feedback control law ()4.4p . and the coordinate 
transformation ()4.2p . we obtain the state feedback control law k{(j){x), u) making S (5g-ISS with respect to the 
input V. The corresponding contraction metric, with respect to states and inputs, for the control system E is 
given by: 

(rG)(x) = e^(x)G(0(.T))e(.T), 

where 

«^ r 1 00 

(x) = 1 

VsG cos{So + xi)x3 VsGsm{So + xi) 

Since the map ((> does not transform the first and the second coordinates and the metric G is only a function 
of the first coordinate, we have G{(j){x)) = G{x). 



(4.5) 



e(x) 



dx 



5. Discussion 

In this paper we developed a backstepping procedure to design controllers enforcing incremental stability. 
Where before we could apply backstepping to construct stabilizing controllers, we can now apply the results 
in this paper to construct incrementally stabilizing controllers. 
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